Connecting Irreversible to Reversible Aggregation: Time and Temperature 
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We report molecular dynamics simulations of a gel-forming mixture of ellipsoidal patchy particles 
with different functionality. We show that in this model, which disfavors the formation of bond- 
loops, elapsed time during irreversible aggregation — leading to the formation of an extended 
network — can be formally correlated with equilibrium temperature in reversible aggregation. We 
also show that it is possible to develop a parameter-free description of the self-assembly kinetics, 
bringing reversible and irreversible aggregation of loopless branched systems to the same level of 
understanding as equilibrium polymerization. 
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Several natural and synthetic materials, as well as bi- 
ological structures, result from the self-assembly of el- 
ementary units into branched aggregates and networks 
[U, [2, S i, H- This self-assembly process is receiv- 
ing considerable attention in two fast-growing fields: 
supramolecular chemistry P, 0, Q and collective be- 
havior of patchy and functionalized particles [1, 0i 
among the most promising building blocks of new ma- 
terials 0, ■ The process of formation of an extended 
three-dimensional network of bonds connecting indepen- 
dent molecules, proteins or colloidal particles, is named 
gelation and the resulting material a gel [S ■ The 
ratio between the bond energy uq and the thermal energy 
UbT (where ks is the Boltzmann constant and T is the 
temperature) can be used to classify aggregation into two 
broad categories: for strong attraction strength {chemi- 
cal case [ll|, [13, E, [13 ) J bond formation is irreversible 
and the number of bonds continuously grows with time. 
In the case of weak attraction strength {physical case 
[l3|), bonds break and reform while the number of bonds 
progressively reaches its equilibrium value. In the latter 
case, the final structure of the system can in principle be 
predicted with equilibrium statistical mechanics meth- 
ods. 

The value of the ratio uo/ksT separates the two 
classes. Conceptually, any model of physical aggregation 
may be turned into a chemical model by studying its 
properties following a quench to ksT <C uq. Similarly, 
applying temperatures comparable to the bond energy 
turns an irreversible aggregation model into a physical 
one. The idea of a close connection between irreversible 
and reversible aggregation is already contained in the 
early mean- field theoretical work of Stockmayer , con- 
sidering the Smoluchowski's kinetic equations solved in 
the limit of absence of closed bonding loops. In Stock- 
mayer's calculations, at any time t during chemical aggre- 
gation, the distribution of clusters of finite size k {Nk{t)) 
is identical to that found following equilibrium statisti- 
cal mechanics prescriptions, i.e., by maximizing the en- 
tropy with the constraint of a fixed number of bonds. 
The Nk are commonly referred to as Flory-Stockmayer 
(FS) distributions. Later on. Van Dongen and Ernst [18l | 



confirmed that the FS distributions are also solutions of 
the Smoluchowski's equations when bond-breaking pro- 
cesses are accounted for. According to these theoret- 
ical works — based on kinetic equations derived in the 
limit of reaction-controlled rates (as opposed to diffusion- 
controlled [3) — a system forming progressively larger 
and larger loopless branched aggregates evolves in time 
via a sequence of states which are identical to the states 
explored in equilibrium at appropriate values of T. The 
equality in the fraction of formed bonds p (the extent of 
reaction in chemical language) provides the connection 
between t during reversible or irreversible aggregation 
and T in equilibrium. Van Dongen and Ernst [l8| also 
provide an analytic expression for the f-dependence of 
p following a sudden change in the external control pa- 
rameters, offering the first soluble example of reversible 
self-assembly of loopless branched structures. Interest- 
ingly, for functionality two (chain assembly) the solution 
coincides with the mean-field analytic expressions later 
on derived for equilibrium polymerization kinetics [20l |. 
Despite their relevance, to our knowledge the Van Don- 
gen and Ernst predictions have never been tested exper- 
imentally or numerically to assess their limit of validity. 

In this Letter — stimulated by the renewed interest 
in reversible and irreversible self-assembly of nano- and 
micro-particles interacting via specific directional bonds 
0, i, i S i — we provide a stringent test of the sug- 
gested connections between chemical and physical aggre- 
gation, and establish the limits of analytic description 
of the kinetics of formation of branched structures. We 
do so by extending a model for chemical gels, recently 
introduced [2l|, to the corresponding physical case by 
turning the attraction strength between bonding sites fi- 
nite. The model, originally inspired by epoxy-amine step- 
polymerization [24I , represents two types of mutually re- 
active molecules A and B as hard homogeneous ellip- 
soids of revolution whose surface is decorated in a prede- 
fined geometry by Ja and Jb identical reactive sites. The 
study of the chemical version of this model [271 showed 
that formation of closed bonding loops in finite size clus- 
ters is disfavored, possibly due to the non-spherical parti- 
cle shape and the location of the reactive sites. Therefore, 
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it offers to us the possibility to carefully check the Stock- 
mayer [ij] and Van Dongen and Ernst [l^ predictions. 
We can assess how closely — in the absence of loop for- 
mation — the chemical gelation process in a system of 
functionalized or patchy units, as well as the reversible 
evolution from a monomeric to an equilibrium bonded 
state, can be envisaged as a progressive sequence of equi- 
librium states, closely connecting t with T. 

We study via event-driven molecular dynamics simu- 
lations a binary mixture composed of Na = 480 penta- 
functional (/^ = 5) ellipsoids of type A and Nb = 1200 
bifunctional {/b = 2) ellipsoids of type B, so that the 
number Na/a of ^-type reactive sites equals the num- 
ber NbJb of _B-type reactive sites. A-particles have 
mass TO, revolution axis a — lOcr and the other two 
axes b = c = 2(t; B-particles have mass 3.4m, and axes 
a — 20(7, 6 = c = 4(7. Size- and mass-ratio are chosen 
to mimic the values of an epoxy-amine system [2^ . The 
interaction potential is the hard ellipsoid potential Ybe 
supplemented by site-site square-well attractive interac- 
tions Vsw (of strength mq, and width (5 = 0.2cr) between 
pairs of particles of different type. Our unit mass is to, 
the unit energy uq- 2^ is measured in units of the potential 
depth (i.e., fc^ = 1), time t in units of cr(TO/uo)^^^- The 
packing fraction is fixed at = 0.3 and T is varied from 
T = 0.3 to r = 0.065. Event-driven (ncwtonian) molec- 
ular dynamics simulations (in the NVT ensemble) are 
performed using a specifically designed code ,23]. Two 
sites, on particles of different type, form a bond if their 
distance is closer than 8. Clusters are defined as groups 
of bonded particles. The particle shape and the location 
of the attractive sites (see inset of Fig. [T|) is such that 
each site is engaged at most in one bond, ensuring an 
unambiguous definition of the extent of reaction p by the 
ratio between the number of bonds in the system and the 
maximum number of possible bonds (i.e., NaJa)- Unless 
otherwise stated, we study the trajectory of the system 
in configuration space, starting from an initial configura- 
tion with no bonds between particles, at fixed T. Dur- 
ing the simulation, bonds form and break continuously 
with time while the system evolves toward the equilib- 
rium state, characterized by an equilibrium value Peq{T) 
of the extent of reaction. To improve the statistics we 
average over 11 different realizations for each studied T . 
Equilibrium results are obtained from the final part of 
the trajectories, after the equilibration transient is over. 

Figure [T] shows Peq (T) , for all T values where equili- 
bration was achieved within the allocated computational 
time (six months at low T). Within the range 0.2 > 
kBT/uQ > 0.08 the system crosses from a monomeric 
state to a significantly bonded state, with more than 
80% of the bonds formed. The sharp sigmoidal shape 
can be perfectly represented by the independent-bond 
mass-action law [23 | 

Peg ^ g/3(«o-TAS) /-[^^ 
(1 -Peg)2 

where uq and AS" describe, respectively, the energy and 



the entropy change associated to the formation of a single 
bond. 

Figures [2lJ a)- (b) show the equilibrium structure of the 
system and the cluster size distribution Nk at three 
different values of T, corresponding to Peq values be- 
low, at and above percolation (numerically located at 
Pc = 0.505 ± 0.007 [lH). To describe the structure of 
the system we report the center-to-center pair distribu- 
tion function gBBif) for the B particles. To quantify 
the T dependence of Nk we also report the number- 
(/c„) and weight-average size (fc^,) of the finite clusters 
(Fig. mjc)). Data shown in Figs. [IKb)-(c) are very well 
described — without any fitting parameter — both below 
and above percolation by the FS predictions [H, [2l| spe- 
cialized to the /a=5-/b=2 case, confirming that bonding 
loops in finite size clusters can be neglected. Therefore, 
the model allows us to check, without fitting parameters, 
if the evolution of the system during equilibration and 
irreversible aggregation does follow a sequence of equi- 
librium steps, and in particular, how t in chemical gela- 
tion can be associated to a corresponding equilibrium T . 
To this aim, we compare in Figs.[2Ka)-(c) the equilibrium 
quantities with the corresponding quantities evaluated at 
selected elapsed times t during the irreversible bonding 
process. The specific t value is chosen in such a way that 
p{t) = peq{T). For all these quantities, at each t, data 
are identical to those obtained in equilibrium at the same 
extent of reaction, demonstrating that the evolution of 
the system structure and connectivity during irreversible 
aggregation does follow a sequence of equilibrium states. 

Also the equilibration process at finite T, i.e., when 
bond-breaking processes are included, proceeds along a 
sequence of equilibrium states described by the FS the- 
ory. Indeed, Fig. [3] shows Nk{t) for several successive 
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FIG. 1: Temperature dependence of the fraction of bonds 
at equilibrium, peq. Symbols: simulation results. The solid 
line represents peq{T) from Eq. [1] A one-parameter best-fit 
gives AS/kB = 8.52. The percolation threshold pc = 0.5 
is indicated by a horizontal line. The inset shows a graphic 
description of the A and B particles. The centers of the small 
spheres locate the bonding sites on the surface of the hard- 
core ellipsoid. 
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FIG. 2: (a) Partial radial distribution function for B parti- 
cles, gBBir); (b) cluster size distribution Nk and (c) number- 
and weight-average size of the finite clusters (fc„ and fc^). In 
all panels, equilibrium results (closed symbols) are compared 
with results at different times t during irreversible aggrega- 
tion (open symbols, from Ref. [2ll ) such that pit) = peq{T). 
In panel (a), curves at T=0.11 and T=0.0826 are shifted by 
-1-0.15 and -1-0.30, respectively. In (b) and (c) solid lines indi- 
cate FS predictions. 



times following a quench, in conditions of reversible ag- 
gregation, from a high-T unassembled state to a low-T 
assembled one. Data are compared — without fit pa- 
rameters — with the corresponding equilibrium quanti- 
ties, accurately modeled by the FS distributions. Similar 
agreement is observed for all the investigated equilibra- 
tion processes. 

Data in Fig. [3] show that the cluster size distribution 
during the kinetics of equilibration evolves via a sequence 
of FS distributions. Fig. 0] shows the associated time 
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FIG. 3; Distribution of finite size clusters A'^^ for different 
times (and fraction of bonds p) during reversible aggregation 
following a T-jump to T = 0.11. At each time during the 
equilibration process, the cluster size distribution is identical 
to the FS distribution at the same fraction of bonds (solid 
lines) 



dependence of the bond probability p{t) , following a T- 
jump starting either from a high-T unbonded configura- 
tion (p(0) = 0) or from a previously equilibrated config- 
uration characterized by branched aggregates (p(0)_^ 0). 
Lines in Fig. |4] are the theoretical predictions [18l | — 
with the Flory post-gel assumption [l^ for values ofp{t) 
above percolation — i.e. 
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breaking 



^szte (1 +Peg)/(1 -Peg), .v.^j^ 

is the rate constant for breakage of a single bond. Since 
Peg is known from Eq. [Tl the entire equilibration dynam- 
ics in this aggregating branching system only depends 
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which fixes the time scale 



on one parameter, jv^j^g 
of the aggregation process. The theoretical expression 
very well represents the numerical data, except for the 
two lowest studied temperatures (T < 0.07), at which 
extensively bonded states are reached (p > 0.8). Such 
disagreement can be traced to a failure of the Flory post- 
gel assumption and/or to a progressive role of the size- 
dependence of the cluster mobility. The overall agree- 
ment between theory and simulation in Fig. |4] is further 
stressed by the T dependence of the single fit-parameter 
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/\/T (see inset), which as expected 
incorporates both an Arrhenius term and the thermal 
velocity component ^ \/T, entering the attempt rate of 
bond breaking. 

The results of this study clearly indicate that the irre- 
versible evolution of a system of patchy or functionalized 
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FIG. 4: Time dependence of the fraction of bonds p{t) dur- 
ing reversible aggregation at different temperatures (r=0.17, 
0.13, 0.12, 0.11, 0.10, 0.095, 0.09, 0.082, 0.07 from the lower 
to the higher equilibrium value p{oo) = Peg), starting from 
an unbonded configuration (p(0) = 0; open symbols) or from 
a generic equilibrated configuration (p(0) 7^ 0; closed sym- 
bols). Lines are solutions of the Smoluchowski's equation 
for reversible aggregation (Eq. [2]), with p^q from Eq. [1] and 
^site''^"^ as the only fit-parameter. The inset shows the T 
dependence of K^^^.^ 

particles, in which bonding loops can be neglected, can 
be put in correspondence with a sequence of equilibrium 
states. For this class of aggregating systems, it is thus 
possible to convert irreversible-aggregation time into an 
effective T and to envisage the evolution of a chemical 
gel as a progressive cooling of the corresponding physical 
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